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Abstract 
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In supergravity models, the evaporation of light Primordial Black Holes 
(PBHs) should be a source of gravitinos. By considering this process, new 
stringent limits are derived on the abundance of small black holes with ini- 
S^ ■ tial masses less than 10 9 g. In minimal supergravity, the subsequent decay 

of evaporated gravitinos into cascades of non-equilibrium particles leads 
to the formation of elements whose abundance is constrained by observa- 
tions. In gauge mediated supersymmetry breaking models, their density is 
required not to overdose the Universe. As a result, cosmological models 
with substantial inhomogeneities on small scales are excluded. 



1 Introduction 

Primordial Black Holes (PBHs) are a very sensitive cosmological probe for physics 
phenomena occurring in the early Universe. They could be formed by many dif- 
ferent mechanisms, e.g., initial density inhomogeneities pQ, a softening of the 
equation of state [3 El El, collapse of cosmic strings [5], a double inflation sce- 
nario jHj, phase transitions jjj, a step in the power spectrum jH], etc. (see [5] for a 
review) . Constraining the abundance of primordial black holes can lead to invalu- 
able information on cosmological processes, particularly as they are probably the 
only viable probe for the power spectrum on very small scales which remain far 
from the Cosmological Microwave Background (CMB) and Large Scale Structures 
(LSS) sensitivity ranges. To date, only PBHs with initial masses between ~ 10 9 g 
and ~ 10 16 g have led to stringent limits (see e.g. [TOll^) from consideration 
of the entropy per baryon, the deuterium destruction, the 4 He destruction and 
the cosmic-rays currently emitted by the Hawking process J2j- The existence 
of light PBHs should lead to important observable constraints, either through 
the direct effects of the evaporated particles (for initial masses between 10 14 g 
and 10 16 g) or through the indirect effects of their interaction with matter and 
radiation in the early Universe (for PBH masses between 10 9 g and 10 14 g). In 
these constraints, the effects taken into account are those related with known par- 
ticles. However, since the evaporation products are created by the gravitational 
field, any quantum with a mass lower than the black hole temperature should 
be emitted, independently of the strength of its interaction. This could provide 
a copious production of superweakly interacting particles that cannot not be in 
equilibrium with the hot plasma of the very early Universe. 

Following [HI EH C3 and fHl 1110 (but in a different framework and using 
more stringent constraints), this article derives new limits on the mass fraction of 
black holes at the time of their formation {(3 = Ppbh/ Ptot) using the production of 
gravitinos during the evaporation process. Depending on whether gravitinos are 
expected to be stable or metastable, the limits are obtained using the requirement 
that they do not overdose the Universe and that the formation of light nuclei by 
the interactions of 4 He nuclei with nonequilibrium flux of D,T, 3 He and 4 He does 
not contradict the observations. This approach is more constraining than the 
usual study of photo-dissociation induced by photons-photinos pairs emitted by 
decaying gravitinos. It opens a new window for the upper limits on j3 below 
10 9 g. The cosmological consequences of our new limits are briefly reviewed in 
the framework of 3 different scenarios : a blue power spectrum, a step in the 
power spectrum and first order phase transitions. 

2 New limits on the PBH density 

Several constraints on the density of PBHs have been derived in different mass 
ranges assuming the evaporation of only standard model particles : for 10 9 g < 



M < 10 13 g the entropy per baryon at nucleosynthesis was used fJJ to obtain 
(3 < (10 9 g/M), for 10 9 g < M < 10 11 g the production of nn pairs at nucleosyn- 
thesis was used (Hj to obtain (3 < 3 x 10" 17 (10 9 g/M) 1 / 2 , for 10 10 g < M < 10 11 g 
deuterium destruction was used (THj to obtain (3 < 3 x 10~ 22 (M/10 10 g) 1 / 2 , 
for 10 11 g < M < 10 13 g spallation of 4 He was used [201 IS] to obtain (3 < 
3 x 10~ 21 (M/10 9 g) 5//2 , for M^5x 10 14 g the gamma-rays and cosmic-rays were 
used PTII22] to obtain (3 < 10 -28 . Slightly more stringent limits where obtained 
in J2H|, leading to (3 < 10~ 20 for masses between 10 9 g and 10 10 g and in [21] . 
leading to (3 < 10 -28 for M = 5 x 10 11 g. Gamma-rays and antiprotons where also 
recently re-analyzed in |2S] and [2E], improving a little the previous estimates. 
Such constraints, related to phenomena occurring after the nucleosynthesis, ap- 
ply only for black holes with initial masses above ~ 10 9 g. Below this value, 
the only limits are the very weak entropy constraint (related with the photon- 
to-baryon ratio) and the quite doubtful relics constraint (assuming stable black 
holes forming at the end of the evaporation mechanism as described, e.g., in [27]). 

To derive a new limit in the initial mass range Mpi < M < 10 11 g, gravitinos 
emitted by black holes are considered in this work. Gravitinos are expected to 
be present in all local supersymmetric models, which are regarded as the more 
natural extensions of the standard model of high energy physics (see, e.g., J2H] for 
an introductory review). In the framework of minimal Supergravity (mSUGRA), 
the gravitino mass is, by construction, expected to lie around the electroweak 
scale, i.e. in the 100 GeV range. In this case, the gravitino is metastable and 
decays after nucleosynthesis, leading to important modifications of the nucle- 
osynthesis paradigm. Instead of using the usual photon-phot ino decay channel, 
this study relies on the more sensitive gluon-gluino channel. Based on |29j . the 
antiprotons produced by the fragmentation of gluons emitted by decaying grav- 
itinos are considered as a source of nonequilibrium light nuclei resulting from 
collisions of those antiprotons on equilibrium nuclei. Then, 6 Li, 7 Li and 7 Be nu- 
clei production by the interactions of the nonequilibrium nuclear flux with 4 He 
equilibrium nuclei is taken into account and compared with data (this approach 
is supported by several recent analysis |3U] which lead to similar results). The 
resulting Monte-Carlo estimates [2H] lead to the following constraint on the con- 



centration of gravitinos: 71,3/2 < '. 1.1 x 10~ 13 m 3 , 9 where 7723/2 is the gravitino 
mass in GeV. This constraint has been successfully used to derive an upper limit 
on the reheating temperature of the order [22]: Tr < 3.8 x 10 6 GeV. The con- 
sequences of this limit on cosmic-rays emitted by PBHs was considered, e.g., in 
|3*T] . In the present approach, we relate this stringent constraint on the gravitino 
abundance to the density of PBHs through the direct gravitino emission. The 
usual Hawking formula J2] is used for the number of particles of type % emitted 
per unit of time t and per unit of energy Q. Introducing the temperature defined 
by T = hc 3 /(167ikGM) pa (10 13 g)/M GeV, taking the relativistic approximation 
for T s , and integrating over time and energy, the total number of quanta of type 



% can be estimated as: 

n tot _ 27 x 10 24 r Tpi dT r x x 2 dx 
1 = Q^asuGRA J Tl T* L/t e x - (-1) 3 

where T is in GeV, Mpi ps 10~ 5 g, x = Q/T, m is the particle mass and asuGRA 
accounts for the number of degrees of freedom through M 2 dM = —ctsuGRAdt 
where M is the black hole mass. Once the PBH temperature is higher than 
the gravitino mass, gravitinos will be emitted with a weight related with their 
number of degrees of freedom. Computing the number of emitted gravitinos as 
a function of the PBH initial mass and matching it with the limit on the grav- 
itino density imposed by nonequilibrium nucleosynthesis of light elements leads 
to an upper limit on the PBH number density. If PBHs are formed during a 
radiation dominated stage, this limit can easily be converted into an upper limit 
on (3 by evaluating the energy density of the radiation at the formation epoch. 
The resulting limit is shown on Fig. 1 and leads to an important improvement 
over previous limits, nearly independently of the gravitino mass in the interest- 
ing range. This opens a new window on the very small scales in the early Universe. 

It is also possible to consider limits arising in Gauge Mediated Susy Break- 
ing (GMSB) models JH2]- Those alternative scenarios, incorporating a natural 
suppression of the rate of flavor- changing neutral-current due to the low energy 
scale, predict the gravitino to be the Lightest Supersymmetric Particle (LSP). 
The LSP is stable if R-parity is conserved. In this case, the limit is obtained 
by requiring f2 3 / 2i o < &m,o, i-e- by requiring that the current gravitino density 
does not exceed the matter density. It can easily be derived from the previous 
method, by taking into account the dilution of gravitinos in the period of PBH 
evaporation and conservation of gravitino to specific entropy ratio, that : 

P< ^ r 



7Vc<-^2/2 I bl 



'3/2 ~M~ [1 



where N 3 / 2 is the total number of gravitinos emitted by a PBH with initial mass 
M, t eq is the end of RD stage and tf = max(tf orm ,t en( i) when a non-trivial 
equation of state for the period of PBH formation is considered, e.g. a dust-like 
phase which ends at t en d [33]. The limit does not imply thermal equilibrium of 
relativistic plasma in the period before PBH evaporation and is valid even for 
low reheating temperatures provided that the equation of state on the preheating 
stage is close to relativistic. With the present matter density fiufi ~ 0.27 [33] 
this leads to the limits shown in Fig. 2 for two extreme cases: 7713/2 = 10~ 5 GeV 
and 777.3/2 = 10 GeV. These results are close to the previous ones and remain very 
competitive in this mass range. Models of gravitino dark matter with ^3/2,0 = 
Qcdai,o, corresponding to the case of equality in the above formula, were recently 
considered in 
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Figure 1: Constraints on the fraction of the Universe going into PBHs (adapted 
from [TUl ITT]). The two curves obtained with gravitinos emission in mSUGRA 
correspond to m 3 / 2 = 100 GeV (lower curve in the high mass range) and 777.3/2 = 
1 TeV (upper curve in the high mass range) 
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Figure 2: Constraints on the fraction of the Universe going into PBHs (adapted 
from [TQJI^). The two curves obtained with gravitinos emission in GMSB cor- 
respond to m 3 / 2 = 10 -5 GeV (lower curve in the high mass range) and 7713/2 = 10 
GeV (upper curve in the high mass range) 



3 Cosmological consequences 

Our new stringent upper limits on the fraction of the Universe in primordial black 
holes can be converted into cosmological constraints on models with significant 
power on small scales. 

The easiest way to illustrate the importance of such limits is to consider a 
blue power spectrum and to derive a related upper value on the spectral index 
n of scalar fluctuations (P(k) oc k n ). It has recently been shown by WMAP [35] 
that the spectrum is nearly of the Harrison-Zel'dovich type, i.e. scale invariant 
with n ~ 1. However this measure was obtained for scales between 10 45 and 10 60 
times larger that those probed by PBHs and it remains very important to probe 
the power available on small scales. The limit on n given in this paper must 
therefore be understood as a way to constrain P(k) at small scales rather than 
a way to measure its derivative at large scales : it is complementary to CMB 
measurements. Using the usual relations between the mass variance at the PBH 
formation time a nit form) and the same quantity today an (to) 



<T H (tfor m ) ~ <TH(t ) [m^Q) 



n — 1 n — 1 



M H (t Q )\ 6 / M H (t eq ) \ ^ 



M H (t 



form, 



where Mjj(t) is the Hubble mass at time t and t eq is the equilibrium time, it is 
possible to set an upper value on (3 which can be expressed as 



,N" 



(T H (tform) 



fj ~ j _g H^ form' 

where 5 m i n « 0.3 is the minimum density contrast required to form a PBH. The 
limit derived in the previous section leads to n < 1.20 in the mSUGRA case 
whereas the usually derived limits range between 1.23 and 1.31 [HniEZj- In the 
GMSB case, it remains at the same level for 1713/2 ~ 10 GeV and is slightly re- 
laxed ton < 1.21 for m 3 / 2 ~ 100 keV. This substantial improvement is due to 
the much more important range of masses probed and where derived for M ~ 1 g. 

In the standard cosmological paradigm of inflation, the primordial power spec- 
trum is expected to be nearly -but not exactly- scale invariant [41] . The sign of 
the running can, in principle, be either positive or negative. It has been recently 
shown that models with a positive running a s , defined as 

Pit) = P( h) (-) 

are very promising in the framework of supergravity inflation (see, e.g., J4"2"]). 
Our analysis strongly limits a positive running, setting the upper bound at a 
tiny value a s < 2 x 10~ 3 . This result is more stringent than the upper limit 



obtained through a combined analysis of Lya forest, SDSS and WMAP data 
[25] . —0.013 < a s < 0.007, as it deals with scales very far from those probed by 
usual cosmological observations. The order of magnitude of the running naturally 
expected in most models -either inflationary ones (see, e.g., |44j ) or alternative 
ones (see, e.g., [35])- being of a few times 10~ 3 our upper bound should help to 
distinguish between different scenarios. 

In the case of an early dust-like stage in the cosmological evolution [21 El E] , 
the PBH formation probability is increased to (3 > 5 13 ' 2 where 5 is the density 
contrast for the considered small scales. The associated limit on n is strengthened 
ton< 1.19. 



Following [36], it is also interesting to consider primordial density perturba- 
tion spectra with both a tilt and a step. Such a feature can arise from underlying 
physical processes J3S] and allows investigation of a wider class of inflaton poten- 
tials. If the amplitude of the step is defined so that the power on small scales 
is p~ 2 times higher than the power on large scales, the maximum allowed value 
for the spectral index can be computed as a function of p. Figure 3 shows those 
limits, which become extremely stringent when p is small enough, for both the 
radiation-dominated and the dust-like cases. The different values of the gravitino 
mass considered in the first section of this paper are all included within the width 
of the lines. 

Another important consequence of our new limits concerns PBH relics dark 
matter. The idea, introduced in J3H], that relics possibly formed at the end of the 
evaporation process could account for the cold dark matter has been extensively 
studied. The amplitude of the power boost required on small scales has been de- 
rived, e.g., in J3U] as a function of the relic mass and of the expected density. The 
main point was that the "step" (or whatever structure in the power spectrum) 
should occur at low masses to avoid the constraints available between 10 9 g and 
10 15 g. The limit on j3 derived in this work closes this dark matter issue except 
within a small window below 10 3 g. 

Finally, the limits also completely exclude the possibility of a copious PBH 
formation process in bubble wall collisions [H)] . This has important consequences 
for the related constraints on first order phase transitions in the early Universe. 

4 Conclusion and prospects 

If local supersymmetry - supergravity - is the correct extension of the standard 
model, the emission of gravitinos from evaporating primordial black holes leads 
to important constraints on their number density in the - so far nearly unexplored 
- low formation mass range. This allows us to exclude cosmological models with 
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Figure 3: Upper limit on the spectral index of the power spectrum as a function 
of the amplitude of the step. 



too much power on small scales. In particular, a blue power spectrum or a posi- 
tive running spectral index are strongly disfavored. Any mechanism which would 
lead to small scale inhomogeneities with a density contrast above ~ 0.3 beyond 
a tiny fraction of the order of 10~ 12 — 10~ 20 for masses between 10~ 5 and 10 10 g 
are excluded by this analysis. This method offers nontrivial links between the 
inhomogeneity of the early Universe and the existence of (meta-) stable particles. 

It should also be noticed that, as in inflationary cosmology, the equilibrium 
is established only after reheating at T < Tr. If Tr is as low as [29J Tr < 3.8 x 
10 6 GeV, superheavy particles with masses m 3> Tr and superweakly interacting 
particles with interaction cross section a <C T l A[ cannot be in equilibrium but 
can be copiously produced through PBH evaporation. This property opens a 
wide range of possible applications for testing particle theory and cosmological 
scenarios. 
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